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Abstract 

The purpose of this short article is to announce, and briefly describe, a 
Maple package, PARTITIONS, that (inter alia) completely automatically dis- 
covers, and then proves, explicit expressions (as sums of quasi-polynomials) 
for p m (n) for any desired m. We do this to demonstrate the power of "rigor- 
ous guessing" as facilitated by the quasi-polynomial ansatz. 

Keywords: integer partitions 
2010 MSC: 05A17, 11P81. 



1. Introduction 

Recall that a partition of a non-negative integer n is a non-increasing 
sequence of positive integers Ai . . . \ m that sum to n. For example the integer 
5 has the following seven partitions: {5,41,32,311,221,2111,11111}. The 
bible on partitions is George Andrews' magnum opus [1]. 
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We denote by p m (n) the number of partitions of n into at most m parts. 
By a classic theorem [1, p. 8, Thm. 1.4], p m (n) also equals the number 
of partitions of n into parts that are at most m. There is an extensive 
literature concerning formulae for p m (n), including contributions by Cay- 
ley, Sylvester, Glaisher, and Gupta. For additional references and historical 
notes, see George Andrews' fascinating article [2, §3] and Gupta's Tables [8, 
pp. i-xxxix]. For an exhaustive history through 1920, see Dickson [4, Ch. 
3]. 

More recently, George Andrews' student, Augustine O. Munagi, devel- 
oped a beautiful theory of so-called g-partial fractions [11], where the de- 
nominators in the decomposition are always expressions of the form (1 — q r ) s , 
rather than powers of cyclotomic polynomials as is the case with the ordi- 
nary partial fraction decomposition. Accordingly, formulae for p m (n) derived 
from the g-partial fraction decomposition of the generating function are most 
naturally expressed in terms of binomial coefficients. 

It is well-known and easy to see that for any m, p m {n) is a sum of quasi- 
polynomials of periods 1,2,3, ... ,m. A quasi-polynomial of period r is a func- 
tion f(n) on the integers such that there exist r polynomials P\{n), P2{n), . . . ,P r 
such that f(n) = Pi(n) if n = i (mod r). We represent such a quasi- 
polynomial as a list [P\{n), . . . , P r (n)]. 

Thus, e.g., we have, for n > 0, 
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Eqs. (l)-(5) were given in 1856 by Cayley [3, p. 132] in a somewhat 
different form. In 1909, Glaisher [6] presented formulae for p m (n) for m = 

I, 2, ... , 10. In 1958, Gupta [8] extended Glaisher's results to the cases m = 

II, 12. In his 2005 Ph.D. thesis [10], Munagi gave formulae for the cases 
m = 1, 2, . . . , 15. Munagi's formulae were derived with the aid of a Maple 
package he developed, and are of a somewhat different character than earlier 
contributions, as they follow from his theory of g-partial fractions [11]. 

2. The PARTITIONS Maple package 

2.1. Overview 

The purpose of this short article is to announce and briefly describe a 
Maple package, PARTITIONS, that completely automatically discovers and 
proves explicit expressions (as sums of quasi-polynomials) for p m (n) for any 
desired m. So far we only bothered to derive the formulae for 1 < m < 70, 
but one can easily go far beyond. 

Not only that, we can, more generally, derive (and prove!), completely 
automatically, expressions, as sums of quasi-polynomials, for the number of 
ways of making change for n cents in a country whose coins have denomina- 
tions of any given set of positive integers. 

Not only that, we can derive (and prove!), completely automatically, 
expressions (as sums of quasi-polynomials) for D k {n), the number of parti- 
tions of n whose Durfee square has size k, for any desired, (numeric) positive 

integer k. (Recall that the size of the Durfee square of a partition Ai X m 

is the largest k such that A& > k.) 

Not only that, we (or rather our computers (and of course yours, if it has 
Maple and is loaded with our package)) can derive asymptotic expressions, 
to any desired order, for both p m (n) and Dk(n). As far we we know the 
formula for Dk(n) is brand-new, and the previous attempts for the asymptotic 
formula for p m {n) by humans G.J. Rieger [14] and E.M. Wright [16] (of 
Hardy-and- Wright fame) only went as far as 0{n~ 2 ) and 0(n -4 ) respectively. 
We go all the way to O(n _10 °)! (and of course can easily go far beyond). 

Not only that, we implement George Andrews' ingenious way [2, sec. 3] 
to convert any quasi-polynomial to a polynomial expression where one is also 
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allowed to use the integer-part function \n\ • This enabled our computers to 
find Andrews- style expressions for p m (n) for 1 < m < 70. 

All these feats (and more!) are achieved by the Maple package PARTITIONS. 

2.2. Using the PARTITIONS package 

In order to use PARTITIONS, you must have Maple™ installed on your 
computer. Then download the file: 

http : //www. math. rutgers . edu/~zeilberg/tokhniot /PARTITIONS and 
save it as PARTITIONS. Then launch Maple, and at the prompt, enter: 

read PARTITIONS: 

and follow the on-line instructions. Let's just highlight the most impor- 
tant procedures. 

AS100(m,n): shows the pre-computed first 100 terms of the asymptotic 
expression, in n, of p m {n) for symbolic m. 

ASD80(k,n): shows the pre-computed first 80 terms of the asymptotic 
expression, in n, of D k {n) for symbolic k. 

BuildDBpmn(n,M) : inputs a symbol n and a positive integer M and out- 
puts a list of size M whose i-th entry is an expression for Pi{n) as a sum of 
i quasi-polynomials 

Discover AS (m, n, k) : discovers the asymptotic expansion to order k of 
Pm{ n ) (the number of partitions of n into at most m parts) for large n and 
fixed, but symbolic, m. 

DiscoverDAS(k,n,r): discovers the asymptotic expansion to order r of 
Dk(n) (the number of partitions of n whose Durfee square has size k) for 
large n and fixed, but symbolic k. 

Durfee (k,n): discovers (rigorously!) the quasi-polynomial expression, in 
n, for D k (n), for any desired positive integer k. It is extremely fast for small 
k, but of course gets slower as k gets larger. 

Durf eePC(k,n) : does the same thing (much faster, of course!) using the 
pre-computed expressions of Durf ee(k,n) ; for k < 40. 

EvalQPS(L,n,nO): evaluates the sum of the quasi-polynomials in the 
variable n given in the list L at n = hq. 

HRR(n,T): evaluates in floating point the sum of the first T terms of the 
Hardy- Ramanuj an- Rademacher formula for p(n), the number of unrestricted 



4 



partitions of n: 



(lyfl (n-£)) 



\ 



p(n) 



1 



ni(s(h,k)-2nh/k) u 

dn 



ttV2 



k>l 0<h<k 

gcd(/i,fc)=l 



I 



where fc) = ^=1 (f - LfJ " 5) (¥ " LfJ - J)is the Dedekind sum. 



Please be warned that for larger n you need to increase Digits. In order 
to get reliable results you may want to use procedure HRRr(n,T,k). 

pmn(m,n): discovers (rigorously!) the quasi-polynomial expression, in n, 
for p m (n), for any desired positive integer to. It is extremely fast for small 
to, but of course gets slower as to gets larger. 

pmnPC(m,n) : does the same thing (much faster, of course!) using the 
pre-computed expressions of pmn(m,n) ; for to < 70. 

pmnAndrews (m,n) : discovers (rigorously!) the Andrews-style expression, 
in n, for p m (n) for any desired positive integer m. Instead of using quasi- 
polynomials explicitly (that some humans find awkward), it uses the integer- 
part function \n\, denoted by trunc(n) in Maple. 

pn(n): the number of partitions of n, p(n), using Euler's recurrence. It 
is useful for checking, since p n (n) = p(n). 

pnSeq(N) : the list of the first iV values of p(n). The output of pnSeq(50000) : 
can be gotten from 

http : //www. math. rutgers . edu/~zeilberg/tokhniot/oPARTITI0NS9 where 
this list of 50000 terms is called pnTable. 

pSn(S,n,K): the more general problem where the parts are drawn from 
the list S of positive integers. It outputs an explicit expression, as a sum 
of quasi-polynomials, for ps{n), the number of integer partitions of n whose 
parts are drawn from the finite list of positive integers S. K is a guessing 
parameter, that should be made higher if the procedure returns FAIL. 

pmnNum(m,nO) : like pmn(m,n) ; but for both numeric to and nO. The 
output is a number. For m < 70 it is extremely fast, since it uses the pre- 
computed values of p m (n) gotten from pmnPC(m,n) ; . For example to get the 
number of integer partitions of a googol (10 100 ) into at most 60 parts, you 
would get, in 0.02 seconds, the 5783-digit integer, by simply typing 

pmnNum(60, 10**100) ; 

One of us (DZ) posed this is a 100-dollar challenge to the users of the 
very useful Mathoverf low forum. This was taken-up, successfully, by user 
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joro[5], whose computer did it correctly in about 2 hours, using PARI. User 
joro generously suggested that instead of sending him a check, DZ should 
donate it in joro's honor, to a charity of DZ's choice, and the latter decided 
on the Wikipedia Foundation. 

Sample input and output can be gotten from the "front" of this arti- 
cle: 

http : //www. math. rutgers . edu/~zeilberg/mamarim/mamarimhtml/pmn.html 

3. Methodology: Rigorous Guessing 

The idea of deriving formulae for p m (n) and ps{n) with the aid of the 
partial fraction decomposition of the generating function dates back at least 
to Cayley [3]. We ask Maple to convert the generating function 



into partial fractions. Then for each piece, Maple finds the first few terms of 
the Maclaurin expansion, and then fits the data with an appropriate quasi- 
polynomial using undetermined coefficients. The output is the list of these 
quasi-polynomials whose sum is the desired expression for p m {n) or ps(n). 
See the source-code for more details. 

Example. Consider the case m = 4. We have Maple calculate that 



4f 4V (i-<?)(i-<? 2 )(i-<? 3 )(i-g 4 ) 

17/72 59/288 1/8 1/24 1/8 1/32 (1 + q)/9 

~ T^q~ + (1 - q) 2 + (1 - qf + (1 - g) 4 + T+q + (1 + q) 2 + 1 + q + q 2 ' 



At this point we could, as Cayley did, expand each term as a series in q, 
collect like terms, and then the coefficient of q n will be a formula for p^n), 
but why bother? From Sylvester [15] and Glaisher [7], we know that 




or in the case of ps(n), where S = {s\, s^, . . . , Sj}, 
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where each Wj(n) is a quasi-polynomial [Pji(n), Pj 2 (n), . . . , Pjj(n)} of period 
j. Further, Wj(n) is of degree [^y^J > an d arises from those terms of (6) with 
denominator a power of the j-th cyclotomic polynomial. Instead, let us allow 
Maple to guess the correct quasi-polynomials: We know a priori that W\ (n) 
is of the form [a + axn + a 2 n 2 + a 3 n 3 ] and let Maple calculate the (beginning 
of the) Maclaurin series for the terms of (6) that contribute to Wi(n): 
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Similarly, for W2(n), which must be of the form 

[ax + a 3 n, a + a 2 n], 



we find 
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Analogous reasoning yields W 3 (n) = [0, — |, |] and W±(n) = [0, — |, 0, |] . 
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4. Conclusion 



The present approach uses very naive guessing to discover, and prove 
(rigorously!), formulas (or as Cayley and Sylvester would say, formulas) for 
the number of partitions of the integer n into at most parts m parts for 
m < 70, and of course, one can easily go far beyond. The core of the idea 
goes back to Arthur Cayley, and is familiar to any second-semester calculus 
student: partial fractions! But dear Arthur could only go so far, so his 
good buddy, James Joseph Sylvester, designated a sophisticated theory of 
"waves" [15] that facilitated hand calculations, which were later dutifully 
carried out by J. W. L. Glaisher in [7]. But, with modern computer algebra 
systems (Maple in our case), one can go much further just using Cayley's 
original ideas. 
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